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J5 (54) Title: METHOD AND APPARATUS FOR NON-INVASIVE ULTRASONIC IMAGING OF HARD TISSUE 

(57) Abstract: A method for non-invasive ultrasonic imaging of hard tissue, comprises the sleps of scanning the biological tissue 
with ultrasonic enei^y transmitted from and received at a plurality of known transmittal and reception locations about the biological 
tissue, to obtain a corresponding plurality of sets of digital data, and processing the sets of digital data to produce an image of a 

2 cross-section of the biological tissue. A corresponding apparatus comprises means for scanning the biological tissue with ulu^onic 
energy from a plurality of transmittal and reception locations about the biological tissue, means for determining the transmittal and 
reception locations, means for generating a plurality of sets of digital data, the sets of digital data representing parameters of the 

^ uluasonic energy at each of the plurality of transmittal and reception locations, and means for processing the plurality of sets of 
digital data to produce a cross-sectional image of the biological tissue. 



wo 02/085178 PCT/IL02/0031 1 

I 

METHOD AND APPARATUS FOR NON-rNVASIVE ULTRASONIC IMAGING OF HARD 

TISSUE 

CROSS-REFERENCE TO RELATED APPLICATIONS 
* 5 This application is entitled to the benefit of U.S. Provisional Application No. 60/284,918 

filed April 20, 2001. 

FIELD AND BACKGROUND OF INVENTION 

The present invention relates, in general, to methods and apparatuses for diagnostic 

10 imaging, and in particular, it discloses a method and apparatus for non-invasive ultrasonic 
imaging of hard tissue. 

Diagnostic imaging of hard tissue has numerous practical uses in various medical fields. 
In the fields of Dentistry, Dental Surgery and Implantology, for example. X-ray and CT imaging 
are extensively used for imaging the human upper and lower jaws. However, such imaging 

1 5 techniques suffer from several significant disadvantages. In order to illustrate the deficiencies of 
current day techniques for imaging hard tissue, let us consider, for example, existing jaw 
imaging techniques used today in dental implant surgery. 

In many mamimals, including humans, the jaws (upper and lower) comprise several layers 
of tissue. FIGURE 1 shows a high-level schematic sectional view of a human lower jaw, or 

20 mandible 30. The external layer of mandible 30 comprises the gum 32, or the mucoperiosteal 
tissue covering the jawbone. Beneath gum 32 is a layer of cortical (compact) bone 34, which is 
normally dense bone tissue. Beneath cortical bone 34 lies an area of trabecular bone 36, which is 
normally bone tissue softer than cortical bone. Within the area of trabecular bone and along the 
mandibular jaw, runs the mandibular canal 38 carrying the inferior alveolar nerve 39. 

25 Mandibular canal 38 is an elongated tubular cavity of varying density, usually comprising dense 
(cortical type bone) borders. In cases where dense borders are not present, mandibular canal 38 is 
a conduit within a sponge-like matrix. In the latter case, mandibular canal cavity 38 may be 
distinguished from the cavities of the trabecular environment 36 by the mandibular canal's 
regular shape, i.e. an elongated tubular cavity. All cavities inside the cross-section of mandible 

30 30 are normally filled with fluids. The upper region of the mandible forms the alveolar ridge 40 
in which teeth 42 are normally situated. 

FIGURE 2a shows a low-level schematic sectional view of mandible 30. As mentioned 
hereinabove, in a normal situation alveolar ridge 40 comprises sockets housing teeth 42. 
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FIGURE 2b shows a sectional view of mandible 30 after the loss of a tooth, for example, due to 
tooth extraction. It is well known that following the loss of a tooth, the residual socket in alveolar 
ridge 40 regenerates and fills with hard tissue. The procedure of replacing a tooth with an 
implant-supported prosthesis has become very common and widely used. The process of fixture 

5 (implant) placement entails exposing the bone (raising the mucoperiosteal flap), and drilling a 
receptive site for the fixture. FIGURE 2c shows a sectional view of mandible 30 with a drill 70 
drilling into the alveolar ridge 40. The common procedure prior to installation of a dental fixture 
involves, inter alia, drilling into the cortical 34 and trabecular 36 bone of mandible 30 using a 
drill 70, in order to prepare a socket in which a fixture will be anchored. If an implantologist 

10 drilling into mandible 30 is not aware of the topology and internal structure thereof (e.g. precise 
location of mandibular canal 38) then during the drilling procedure drill 70 may contact and 
damage mandibular canal 38 and inferior alveolar nerve 39 situated therein, or drill outside the 
bone boundaries. Such damage can ultimately lead to severe pain, hemorrhage and even local 
paralysis, among other undesired consequences. For this reason, it is widely acknowledged that 

1 5 an implantologist must obtain an internal image of the mandible prior to performing the implant 
surgery. 

Similarly, when planning implant surgery on a human upper jaw, or maxilla, an 
implantologist will also need to obtain an internal image of the jaw. FIGURE 3a shows a 
. schematic sectional view of a maxilla 50. The maxilla comprises ah external layer of gum 52, 

20 under which lies a layer of cortical bone 54, enclosing an area of trabecular bone 56. Under 
trabecular bone 56 lies another layer of cortical bone 58 which serves as the floor of the sinal 
and\or nasal cavities 60. In a normal (healthy) maxilla, a tooth 62 is anchored in the cortical 54 
and trabecular 56 bone. FIGURE 3b shows a sectional view of maxilla 50 after the loss of tooth 
62, for example, due to tooth extraction. When preparing a receptive site for a fixture in maxilla 

25 50, the implantologist uses a drill 72 in order to drill a hole through the cortical 54 and trabecular 
56 bone. However, if the implantologist is not aware of the internal structure and shape of the 
jawbone (e.g. the location of sinal and\or nasal cavities 60), the implantologist may perforate and 
damage cortical floor 58. Such perforation may lead to a serious sinal infection and hemorrhage, 
as well as other undesired consequences. It is therefore common practice to obtain an internal 

30 image of the maxilla prior to the implant surgery. 

The most common technique currently used in Implantology for imaging the lower andXor 
upper jaw is Panoramic X-ray Radiography. FIGURE 4 is an example of a panoramic x-ray 
image of human upper and lower jaws. In the lower jaw, mandibular canal 38 (black) appears 
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inside trabecular bone area 36 (gray). A cortical bone layer 34 (white) encloses trabecular bone 
area 36. In the upper jaw, sinai and nasal cavities 60 (black) appear above cortical floor 58 
(white), which is above trabecular bone area 56 (gray), under which lies cortical bone layer 54 
(white). 

5 The panoramic X-ray technique suffers from some significant shortcomings. First, it is 

well established that X-ray radiation is hazardous to the health of the patient. Second, panoramic 
X-ray produces a two-dimensional image of the jaw, which is perpendicular to a cross-section of 
the jaw. This limitation makes the panoramic image unreliable for guiding the implantologist to 
drill within the bone boundaries and within a safe distance from the mandibular canal, or the 

10 sinalXnasal cavities. The panoramic image is inherently distorted and inaccurate because it 
projects the three-dimensional jaw into a two-dimensional image. This image is therefore 
unreliable also for assessing the depth of the bone tissue available for drilling and preparing a 
fixture. Third, the image is not taken chair-side and consequently, panoramic X-ray does not 
allow for real-time monitoring of implant procedures. All of the above disadvantages make the 

1 5 panoramic X-ray image a hazardous, imprecise, and unreliable imaging solution. 

Another imaging technique used in Implantology, though less common, is Computerized 
Tomography (CT). FIGURE 5a is an example of a sectional CT image of a toothless mandible 
30. Cortical layer 34 (white) encloses trabecular area 36 (gray), which surrounds mandibular 
canal 38 (black ellipse). FIGURE 5b is an example of a toothless maxilla 50. Cortical layer 54 

20 (white) covers trabecular area 56 (gray). Cortical floor 58 (while) borders sinal and nasal cavities 
60 (black). 

A CT image of an upper or lower jaw provides a sectional view of the jaw, and is less 
distorted than panoramic radiography. However, CT involves a substantially higher dosage of 
X-ray radiation than conventional radiography, and therefore poses a significantly greater risk to 
25 the health of the patient. Furthermore, CT equipment is very expensive and is only rarely found 
inside the clinic of the implantologist. CT can definitely not provide a chair-side imaging 
solution. 

The popularity of ultrasonic medical diagnostic systems has significantly risen in recent 
years. In contrast to X-ray and CT systems, ultrasonic systems have the advantage of not 
30 exposing the patient or doctor to hazardous ionizing radiation, and are generally more compact 
and economical. In the field of Dentistry and Dental Implantology various ultrasonic diagnostic 
and measurement systems are known. 
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International patent application PCT/ILOO/00341 publication no. WO 01/00102 entitled 
"Alveolar Bone Measurement System" (hereunder "ASMS"), which is ftilly incorporated herein 
by reference, discloses an ultrasound system for assessment of distance between an area of 
interest and a known location of a non-bone canal for use in drilling an implant receiving cavity 
5 in the alveolar bone of a human subject's posterior mandible or posterior maxilla. ABMS 
comprises an ultrasound probe capable of being introduced at the area of interest and 
transceiving pulse echo ultrasound signal to the alveolar bone and therefrom and an electronic 
circuitry for processing the ultrasound signal and providing an indication of the remaining 
alveolar bone distance between the ultrasound probe and a canal within the alveolar bone. 

10 However, ABMS still comes short of answering the needs of the dental implantologist for 

the following reasons. First, measurement of time-of-flight (TOF) from a location on the surface 
of the alveolar bone to a non-bone canal inside the jaw is, in reaHly, impracticable or at least very 
imprecise due to the high level of attenuation and scattering inside the jawbone. Although the 
application further discloses an improved method in which a second TOF measurement is taken 

15 after drilling a bore of known depth, the improvement is still subject to the aforementioned 
attenuation and scattering problem, and moreover, as mentioned in the ABMS patent application 
itself (page 2 line 1) it is of significance that the condition of the jaw be assessed prior to drilling. 
Second, ABMS relies on the so-called "average velocity of ultrasound within bone tissue, as 
known per se" (page 4, line 14). However, it is known that the velocity of ultrasound may vary 

20 from patient to patient, and from bone to bone within a certain patient, and even in different 
regions of a certain bone. Thus, even in the case that ABMS manages to take a precise TOF 
measurement, it will still not be able to calculate the precise distance from the probe to the canal 
due to an error in the velocity of ultrasound. Third, ABMS does not disclose any mechanism or 
procedure for ensuring that the echo which is supposedly from the canal and on which the 

25 distance measurement is based, is really from the canal and not from another reflector inside the 
jawbone. Fourth, ABMS is limited to measuring the distance from the canal to the alveolar bone, 
and does not teach how to measure the distance between the canal and the buccal and lingual 
walls of the jawbone, which is of significant importance to the implantologist, e.g. in order to 
determine an optimal angle of drilling. Lastly, being a measurement system rather than an 

30 imaging system, the most ABMS can provide is a numerical distance measurement from the 
location where the probe is located to a canal within the bone, but no implantologist will suffice 
with a mere numerical value as a basis for planning or performing a drill into a jaw. 
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German Patent No. DE 19921279 (hereunder "the '279 patent"), which is fully 
incorporated herein by reference, discloses a surgical instrument for drilling into a bone, the 
instrument comprising an ultrasonic transducer for transmitting and receiving ultrasonic waves. 
The transducer is connected to a device which generates signals according to the intensity and 
5 TOF of ultrasonic energy received by the transducer, and these signals provide measurements for 
determining the characteristics of the bone in the direction of transmission. The '279 patent 
suffers from limitations similar to those of ABMS, namely, impracticable or imprecise 
measurement and insufficient information to the implantologist. 

U.S. Registered Patent No. 6,030,221 entitled "Ultrasonic Apparatus and for Precisely 

10 Locating Cavitations within Jawbones and the Like" (hereunder "the *221 patent")* which is 
fiilly incorporated herein by reference, discloses an apparatus which generates an ultrasonic 
pulse and passes the pulse through the jawbone of a human. The pulse is detected by an 
ultrasonic receiving unit. Attenuations in the amplitude of the pulse are detected and displayed 
on a color monitor. The color monitor allows the detection of cavitations by interpreting color 

15 codes in a 4x4 matrix displayed on the monitor. 

U.S. Registered Patent No. 6,086,538 entitled "Methods and Apparatus for Evaluation of 
Bone Condition" (hereunder *the '538 patent") discloses a method of evaluating the status of 
bone tissue, useful in the diagnosis of osteoporosis, in which a calcaneus is scanned in 
through-transmission mode, and a characteristic of ultrasound, such as the speed-of-sound or 

20 attenuation, is measured in different locations. The location of a circular (as seen from the side) 
area of reduced attenuation inside the calcaneus is derived from the ultrasound measurements, 
and finally the status of the examined bone tissue is evaluated based on the measurements which 
were taken in that circular area. 

Both the '221 patent and the '538 patent concentrate on the problem of assessing the 

25 quality or health of the bone under examination rather than providing an image of the bone for 
guidance in a medical procedure. As a result, both these patents provide only a two-dimensional 
attenuation map of the examined bone which is perpendicular to a cross-section of the bone, and 
furthermore, include no mechanism or procedure for precisely calculating the distance of the 
detected cavitations (in the case of the '221 patent) or the circular area of reduced attenuation (in 

30 the case of the '538 patent) in relation to a reference point of interest. As mentioned hereinabove 
in connection to the panoramic X-ray technique, a lateral two-dimensional image is unreliable 
for guiding the implantologist to drill within the bone boundaries and within a safe distance from 
the mandibular canal, or the sinal/nasal cavities. 
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German Patent No. DE 4205360 (herevinder "the *360 patent"), which is fully 
incorporated herein by reference, discloses an ultrasonic measuring gauge for determining 
jawbone width. U.S. Registered Patent No. 5,427,105 entitled "Measuring Procedure for the 
Thickness of the Mucous Membrane of an Alveolar Process" (hereunder "the M05 patent"), 
5 which is fully incorporated herein by reference, discloses an ultrasonic method for measuring the 
thickness of the mucous membrane in tlie region of the jawbone ridge. Neither the '360 patent 
nor the M05 patent comprise any mechanism for scanning the bone being examined. Neither 
patent provides an image of the bone, nor any information regarding the internal structure of the 
examined jawbone. 

10 U.S. Registered Patent No. 5,564,423 entitled "Ultrasonic Measurement System for the 

Determination of Bone Density and Structure" (hereunder "the *423 patent"), which is fully 
incorporated herein by reference, discloses an electronic system for measuring the density and 
structure of bone, equipped with ultrasonic calipers designed to be applied to a segment of the 
human body (for example, a finger) containing bone tissue to be examined. The ultrasonic 

15 calipers include a transmitting transducer and a receiving transducer, which enable measuring 
TOP in the bone tissue based on a through-transmission method. The system provides an 
indication of the density and structure of the bone tissue based on the measured TOP. Since the 
'423 patent relies on TOP measurement, it suffers from the same impracticability and inaccuracy 
problems mentioned above in connection to ABMS. Likewise, the '423 patent also does not 

20 comprise any mechanism for scanning the examined bone tissue, and does not provide the 
location and image of internal structures within the bone tissue. 

Thus, none of the above solutions provides an economic, radiation free, real-time, 
chair-side imaging tool to the implantologist. Ongoing monitoring of the drilling process allows 
for depth and angulation corrections on the fly. There is thus a widely recognized need for, and it 

25 would be highly advantageous to have a method and apparatus for non-invasive imaging of hard 
tissue, devoid of the above limitations. 

SUMMARY OF THE INVENTION 

The present invention discloses a method and apparatus for non-invasive ultrasonic 
30 imaging of hard tissue. 

According to the present invention, there is provided in a first embodiment, a method of 
ultrasonic imaging of a biological tissue, comprising the steps of scanning the biological tissue 
with ultrasonic energy transmitted from and received at a plurality of known transmittal and 
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reception locations about the biological tissue, to obtain a corresponding plurality of sets of 
digital data, and processing the sets of digital data to produce an image of a cross-section of the 
biological tissue. 

According to yet another feature of the first embodiment of the method of the present 
5 invention, the step of scanning includes applying at least one scanning mode selected from the 
group consisting of through-transmission mode and pulse-echo mode. 

According to yet another feature of the first embodiment of the method of the present 
invention, the step of scanning includes using al least one scanning mechanism selected from the 
group consisting of mechanical scanning and electronic scanning. 
10 According to yet another feature of the first embodiment of the method of the present 

invention, the step of scanning includes the substeps of disposing at least one ultrasonic 
transducer in proximity to the biological tissue, determining a location of the at least one 
ultrasonic transducer, transmitting ultrasonic energy into the biological tissue using the at least 
one ultrasonic transducer, the ultrasonic energy being reflected from the biological tissue, and 
15 receiving the ultrasonic energy reflected from the biological tissue using the at least one 
ultrasonic transducer. 

According to additional features in the first embodiment of the method of the present 
invention, the method ftirther includes changing the location of the at least one ultrasonic 
transducer, and repeating the steps of: determining a location of the at least one ultrasonic 

20 transducer, transmitting ultrasonic energy into the biological tissue using the at least one 
ultrasonic transducer, the ultrasonic energy being reflected from the biological tissue, and 
receiving the ultrasonic energy reflected from the biological tissue using the at least one 
ultrasonic transducer for a desired number of repetitions. 

According to a second embodiment of the method of the present invention, the step of 

25 scanning includes the substeps of disposing a plurality of ultrasonic transducers in proximity to 
the biological tissue, determining a location of each of the plurality of ultrasonic transducers, 
transmitting ultrasonic energy into the biological tissue, using at least one first ultrasonic 
transducer selected from the plurality of ultrasonic transducers, the ultrasonic energy propagating 
through the biological tissue, and receiving the through propagating ultrasonic energy, using at 

30 least one second ultrasonic transducer selected from the plurality of ultrasonic transducers, the at 
least one second ultrasonic transducer being different from the at least one first ultrasonic 
transducer. 



wo 02/085178 



PCT/IL02/00311 



8 

According to additional features in the second embodiment of the method of the present 
invention, the method further includes the substeps of varying the selection of at least one 
ultrasonic transducer from the group consisting of the at least one first ultrasonic transducer and 
the at least one second ultrasonic transducer, and repeating the steps of: determining a location of 
5 each of the plurality of ultrasonic transducers, transmitting ultrasonic energy into the biological 
tissue, using at least one first ultrasonic transducer selected from the plurality of ultrasonic 
transducers, the ultrasonic energy propagating through the biological tissue, and receiving the 
through propagating ultrasonic energy, using at least one second ultrasonic transducer selected 
from the plurality of ultrasonic transducers, the at least one second ultrasonic transducer being 

10 different from the at least one first ultrasonic transducer, for a desired number of repetitions. 

According to a third embodiment of the method of the present invention, the step of 
scanning includes the substeps of: (a) disposing a plurality of ultrasonic transducers in proximity 
to the biological tissue, (b) determining a location of each of the plurality of ultrasonic 
transducers, (c) transmitting ultrasonic energy into the biological tissue, using at least one 

15 ultrasonic transducer selected from the plurality of ultrasonic transducers, the ultrasonic energy 
being reflected fi-om the biological tissue, (d) receiving the ultrasonic energy reflected from the 
biological tissue using the same at least one selected ultrasonic transducer used for transmitting, 
(e) varying the selection of the at least one selected ultrasonic transducer used for transmitting 
and receiving the ultrasonic energy; and repeating steps (c) through (e) for a desired number of 

20 repetitions. 

According to the present invention, the substep of generating a spectral function further 
includes the substep of representing each spectral function by a first representative value, all 
such first representative values forming a first array of representative values. 

According to another feature in the method of the present invention, the first 
25 representative value is selected from the group consisting of a minimum, a maximum, an 
average, a root mean square (RMS), and a total sum of amplitudes. 

According to yet another feature in the method of the present invention, the substep of 
representing each spectral function by a first representative value further includes the substeps of 
identifying at least one distinct local variability in the first array of representative values, and 
30 calculating the value of a dimension related to an internal structure within the biological tissue, 
based on the at least one distinct local variability. 
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According to yet another feature in the method of the present invention, the dimension 
related to an internal structure within the biological tissue is selected from the group consisting 
of a width, a length, a depth, a height, a diameter and a thickness. 

According to yet another feature in the method of the present invention, the internal 
5 structure includes a cavity. 

According to yet another feature in the method of the present invention, the step of 
generating a spectral function further. includes the substeps of representing each spectral function 
by a second representative value, all such second representative values forming a second array of 
representative value, and producing a synthetic array of representative values, based on the first 
10 and second arrays of representative values. 

According to yet another feature in the method of the present invention, the step of 
generating a spectral function further includes the substeps of identifying at least one distinct 
local variability in the synthetic array of representative values, and calculating the value of a 
dimension related to an internal structure within the biological tissue, based on the at least one 
1 5 distinct local variability. 

According to yet another feature in the method of the present invention, the substep of 
representing each spectral function by a first representative value further includes the substeps of 
identifying at least one extreme value in the first array of representative values, measuring a first 
attenuation coefficient inside the biological tissue, measuring a second attenuation coefficient 
20 inside an internal structure within the biological tissue, measuring the value of a first dimension 
related to the biological tissue, andcalculating the value of a second dimension related to the 
internal structure, based on the at least one identified extreme value, the first attenuation 
coefficient, the second attenuation coefficient, and the value of the first dimension. 

According to the present invention, in a fourth embodiment of the method, the step of 
25 processing includes the substeps of deriving, from each of the sets of digital data, a property 
value related to an ultrasonic property of the biological tissue, measuring, for each of the sets of 
digital data, a dimension value of a dimension related to the biological tissue and to same set of 
digital data, and calculating, for each of the sets of digital data, a set ratio between the property 
value and the dimension value. 
30 According to the present invention, in the fourth embodiment of the method, the step of 

processing further includes the substeps of forming an array of ratios based on the set ratios 
calculated for all the sets of digital data, identifying at least one distinct local variability in the 
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array of ratios, and calculating the value of a dimension related to an internal structure within the 
biological tissue, based on the at least one identified distinct local variability. 

According to the present invention, in the various embodiments of the method, the step 
of processing includes deriving, from each of the sets of digital data, at least one value related to 
5 an ultrasonic property of the biological tissue. 

According to the present invention, in the various embodiments of the method the 
ultrasonic property includes ultrasonic attenuation. 

According to the present invention, in the various embodiments of the method, the step of 
processing includes the substep of generating, for each of the sets of digital data, a spectral 
10 function of the corresponding to the received ultrasonic energy. 

According to the present invention, in the various embodinients of the method, the 
substep of generating a spectral function includes applying a feist Fourier transform (FFT) 
algorithm. 

According to the present invention, in the various embodiments of the method, the step of 
1 5 processing includes applying a Radon transform algorithm. 

According to one feature of all embodiments of the method of the present invention, the 
biological tissue includes a jaw. 

According to another feature of all embodiments of the method of the present invention, 
the biological tissue includes a bone. 
20 According to yet another feature of all embodiments of the method of the present 

invention, the bone is selected from the group consisting of a mandible and a maxilla. 

According to the present invention, there is provided an apparatus for ultrasonic imaging 
of a biological tissue, comprising means for scanning the biological tissue with ultrasonic energy 
from a plurality of transmittal and reception locations about the biological tissue, means for 
25 determining the transmittal and reception locations, means for generating a plurality of sets of 
digital data, the sets of digital data representing parameters of the ultrasonic energy at each of the 
plurality of transmittal and reception locations, and means for processing the plurality of sets of 
digital data to produce a cross-sectional image of the biological tissue. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to the 
accompanying drawings. With specific reference now to the drawings in detail, it is stressed that 
the particulars shown are by way of example, and for purposes of illustrative discussion of the 
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preferred embodiments of the present invention only, and are presented in the cause of providing 
what is believed to be the most useful and readily understood description of the principles and 
conceptual aspects of the invention. In this regard, no attempt is made to show structural details 
of the invention in more detail than is necessary for a fundamental understanding of the 
5 invention, the description taken with the drawings making apparent to those skilled in the art 
how the several forms of the invention may be embodied in practice. 

In the drawings: ^ 

FIGURE 1 shows a high-level schematic sectional view of a mandible; 

FIGURE 2a shows a low-level schematic sectional view of a mandible; 
10 FIGURE 2b shows a sectional view of a mandible after the loss of a tooth; 

FIGURE 2c shows a sectional view of a mandible with a drill drilling into the alveolar 

ridge; 

FIGURE 3a shows a schematic sectional view a maxilla; 

FIGURE 3b shows a sectional view of a maxilla after the loss of a tooth; 
1 5 FIGURE 4 is an example of a panoramic x-ray image of a mandible and a maxilla; 

FIGURE 5a is an example of a sectional CT image of a toothless mandible; 

FIGURE 5b is an example of a sectional CT image of a toothless maxilla; 

FIGURE 6 illustrates a preferred embodiment apparatus of the present invention; 

FIGURE 7 illustrates the preferred method of operation of the preferred embodiment 
20 apparatus of the present invention; 

FIGURES 8a - 8c illustrate an example of a vertical, parallel scanning movement; 

FIGURES 9a - 9c illustrate an example of a horizontal, parallel scanning movement; 

FIGURE 10 illustrates an example of three spectral functions of an ultrasonic signal 
transmitted through a jaw; 
25 FIGURE 1 1 illustrates an alternative embodiment apparatus of the present invention; 

FIGURE 12 illustrates the preferred method of operation of the alternative embodiment 
apparatus of the present invention; 

FIGURES 13a- 13c illustrate an example of a vertical scanning movement; 

FIGURES 14a - 14b illustrate an example of a horizontal scanning movement; 
30 FIGURE 15 illustrates an example of two travel paths of ultrasonic pulses propagating 

through a maxilla; 

FIGURE 16 illustrates yet another alternative embodiment apparatus of the present 
invention; 
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FIGURE 17 illustrates still another alternative embodiment apparatus of the present 
invention 

FIGURE 18 illustrates the preferred method of operation of another alternative 
embodiment apparatus of the present invention; 

5 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention discloses a method and apparatus for non-invasive ultrasonic 
imaging of hard tissue. 

Referring now to FIGURE 6, a preferred embodiment of the present invention, 

10 designated 100, v^ill be described. Apparatus 100 is an apparatus for non-invasive ultrasonic 
imaging of hard tissue. Apparatus 100 comprises a controller 102 coupled to a pulser\receiver 
104 and to a location monitoring device 116. Pulser\receiver 104 is coupled to a first transducer 
106, to a second transducer 108 and to an analog-to-digital converter (ADC) 110. ADC 110 is 
coupled to a computer 112, which is coupled to a display device 114. Location monitoring 

15 device 116 is coupled to first transducer 106, to second transducer 108, and to computer 112. 

Referring now to FIGURE 7, the method of operation of apparatus 100 will be described. 
Apparatus 100 operates according to a through-transmission method, as follows. A user couples 
first transducer 106 to an initial location on a surface of an organ being examined, for example, 
on a first surface 32 of a mandible 30. First surface 32 can be, for example, a surface inside the 

20 mouth cavity, such as the buccal or lingual gum surface of the jaw, or a surface outside the 
mouth cavity, such as the surface of the cheek or the chin. Alternatively, in case the jawbone is 
exposed (e.g. by raising the mucoperiosteal flap) first transducer 106 can be coupled directly to a 
bone surface of the jaw. In FIGURE 7, for example, first surface 32 is a toothless buccal gum 
surface, although the same procedure can be performed on a surface under which a tooth is 

25 present. 

The user couples second transducer 108 to an initial location on a surface opposite to 
transducer 106. In FIGURE 7, for example, second transducer 108 is coupled to a second surface 
32* of mandible 30, which is opposite to first surface 32 to which first transducer 106 is coupled. 
If, for example, first surface 32 is the buccal gum surface of mandible 30, then second surface 
30 32* can be the lingual surface thereof; if first surface 32 is the lingual surface of mandible 30, 
then second surface 32' can be, for example, the buccal gum surface or cheek surface opposite 
thereto. Like first transducer 106, second transducer 108 can be coupled directly to a bone 
surface of the jaw, in case the jawbone is exposed. 
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In order to improve receipt by second transducer 108 of ultrasonic signals emitted by first 
transducer 106, the user preferably aligns the two transducers so that they face one another on a 
predetermined common axis CXio6. In order to further improve penetration of ultrasound signals 
into mandible 30, the user applies a coupling material 74 between first transducer 106 and first 
5 surface 32, and between second transducer 108 and second surface 32'. In some cases, the fluids 
naturally present in the mouth of the patient (e.g. saliva) will be sufficient to serve as coupling 
material 74. In other cases, the user can use another material for coupling, for example, a 
non-toxic ultrasound coupling gel. 

When the transducers are properly coupled, the user starts the operation of apparatus 100 

10 using controller 102, The controller commands location monitoring device 116 to determine the 
initial location of first transducer 106 and the initial location of second transducer 108. Device 
116 determines the initial locations using a conventional measuring method and sends the results 
to computer 112. The controller further commands pulser\receiver 104 to generate electrical 
pulses and to send them to first transducer 106. First transducer 106 converts the electrical pulses 

15 to ultrasonic pulses and emits them toward first surface 32. The emitted ultrasonic pulses 
partially penetrate first surface 32 via coupling material 74. Some of the penetrating ultrasonic 
pulses then travel inside mandible 30 through cortical bone layer 34, trabecular bone 36, opposite 
cortical bone layer 34' and second (opposite) surface 32*, eventually reaching second transducer 
108 via coupling material 74. Second transducer 108 receives ultrasonic pulses which propagated 

20 through mandible 30, converts them into analog electrical pulses and sends the electrical pulses 
to ADC 110 via pulser\receiver 104. The ADC samples the electrical pulses which represent the 
received pulses and sends these samples to computer 112. The computer records, as a first set of 
results, the samples which represent the received pulses along with the initial location of the 
emitting transducer and the initial location of the receiving transducer. 

25 The user scans mandible 30, for example, by moving first transducer 106 along first 

surface 32 whilst also moving second transducer 108 along second surface 32' in parallel to the 
first transducer. FIGURES 8a - Sc illustrate an example of a vertical, parallel scanning 
movement. FIGURE 8a shows the initial locations of transducers 106 and 108 which are 
relatively far from mandibular canal 38. FIGURE 8b shows the locations of the two transducers 

30 after having been moved, in parallel, approximately halfway closer to mandibular canal 38. 
FIGURE 8c shows the locations of the two transducers after having been moved further in 
parallel, so that mandibular canal 38 lies on the ultrasonic travel path (depicted as a broken-line 
arrow) between the two transducers. FIGURES 9a - 9b illustrate an example of a horizontal, 
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parallel scanning movement. First transducer 106 is coupled to alveolar ridge 40 of mandible 30 
and second transducer 108 is coupled to the floor of the mandible, opposite to and facing first 
transducer 106. FIGURE 9a shows the initial locations of the transducers. The transducers are 
initially placed so that the travel path (depicted as a broken-line arrow) of ultrasonic pulses does 
5 not pass through mandibular canal 38. FIGURE 9b shows the location of the two transducers 
after having been moved, in parallel, toward the middle of the jaw, so that mandibular canal 38 
lies on the path between first transducer 106 and second transducer 108. 

Alternatively, instead of a parallel scanning motion (FIGURES 8 and 9), the user may 
also scan mandible 30 by moving only one transducer whilst the other transducer remains 

10 stationary, by moving the two transducers in opposite directions or any other desired scanning 
pattern. Apparatus 100 can further be equipped with automatic mechanical scanning means, for 
example, a slidable mount to which one or both of the transducers are attached and which is 
programmed to execute a desired scanning pattern. Alternatively, apparatus 100 can be equipped 
with electronic scanning means, for example, an array of transducers instead or in addition to 

1 5 transducer 106 andXor transducer 108. 

During the scanning process, ADC 110 samples the ultrasonic pulses received at each 
new location of the transducers, and computer 112 keeps recording sets of samples along with 
corresponding locations of the transducers. Computer 112 is . further programmed to determine 
certain physical characteristics of the hard tissue being examined, based on an analysis of the 

20 recorded sets. 

When ultrasonic pulses travel inside hard tissue, they often pass through layers, areas 
and\or internal structures which possess different acoustic properties. The difference between 
acoustic properties of an internal structure within hard tissue on one hand, and acoustic 
properties of the proximate surroundings of the internal structure on the other hand, makes it 

25 possible for apparatus 100 to detect, locate, measure and image such an internal structure within 
hard tissue. For example, when the hard tissue being examined by apparatus 100 is a mandible, 
penetrating ultrasonic pulses will sometimes propagate through mandibular canal 38 on their 
path from first transducer 106 to second transducer 108. This will be the case, for example, if 
mandibular canal 38 lies on the path between the two transducers (as illustrated in FIGURE 8c). 

30 According to the teachings of the present invention, mandibular canal 38 possesses different 
acoustic properties than trabecular 36 and cortical bone 34. For example, since the mandibular 
canal is typically filled with fluids, it will sometimes be characterized by lower ultrasonic 
attenuation than cortical and trabecular bone. In other cases, the geometry of the mandibular 
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canal may actually cause penetrating ultrasonic pulses to scatter, eventually causing higher 
attenuation than cortical and trabecular bone. The present invention therefore teaches that when, 
during the scanning process, mandibular canal 38 is situated on the path between first transducer 
106 and second transducer 108 (e.g. FIGURES 8c and 9b), consequently computer 112 will be 
3 able to detect a local change (increase or decrease) in amplitude of the received and sampled 
pulses. 

Hence, it is a particular feature of apparatus 100 that it detects and locates an internal 
structure within hard tissue (e.g. the mandibular canal in a mandible), by detecting a local change 
in the amplitude of penetrating ultrasonic pulses. When computer 112 detects such a local 
10 change in the amplitude, it will notify the user (e.g. by issuing an indication on display 114) that 
an internal structure lies on the ultrasonic travel path between first transducer 106 and second 
transducer 108. 

It is another particular feature of apparatus 100 that it can mesisure the depth of an 
internal structure within hard tissue, such as the depth of the mandibular canal in a mandible, 

1 5 based on amplitude measurements. By "depth" is meant the distance from a surface of the organ 
being examined to the internal structure of interest. For purposes of illustration hereunder, 
horizontal depth means the horizontal distance from surface 32 to the nearest border of 
mandibular canal 38, and vertical depth Ly means the vertical distance from alveolar ridge 40 to 
the nearest border of mandibular canal 38 (see FIGURES 8c and 9c). The method of operation of 

20 apparatus 100 in measuring the depth of an internal structure within the hard tissue, for example, 
the depth of mandibular canal 38 within mandible 30, is as follows. The transducers scan 
mandible 30 according to the method described heretofore. As soon as computer 112 first detects 
a characteristic change in the amplitude of penetrating ultrasonic pulses, it means that the pulses 
have first met mandibular canal 38 on their path from first transducer 106 to second transducer 

25 108. Since the computer receives from location monitoring device 116 the location of the 
transducers at any given moment, the computer can therefore calculate the depth of mandibular 
canal as follows. In case of vertical scanning (FIGURES 8a - 8c), vertical depth Ly will be 
calculated as the vertical distance from alveolar ridge 40 to the location of the transducers which 
first introduced a characteristic amplitude change. Similarly, in case of horizontal scanning 

30 (FIGURES 9a - 9c), horizontal depth L, will be calculated as the horizontal distance from 
surface 32 to the location of the first characteristic amplitude change. The measured depth will 
then be displayed on display 114. 
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It is yet another particular feature of apparatus 100 that it can measure the diameter of an 
internal structure of interest within hard tissue, for example, the diameter of mandibular canal 38 
within mandible 30, based on amplitude measurements. The borders of mandibular canal 38 are 
usually elliptic rather than round (see FIGURE 5a), and therefore the diameter may differ 
5 depending on the direction from which it is measured. For purposes of illustration hereunder, 
horizontal diameter Dx means the diameter which lies approximately on the horizontal axis of 
the jaw, and vertical diameter Dy means the diameter which lies approximately on the vertical 
axis thereof (see FIGURES 8c and 9c). The method of operation of apparatus 100 in measuring 
the diameter of an internal structure within hard tissue, for example, the diameter of mandibular 

10 canal 38 in mandible 30, is as follows. The transducers scan mandible 30 according to the 
rnethod described heretofore. When measuring horizontal diameter Dx horizontal scanning 
(FIGURE 9) will be performed, and when measuring vertical diameter Dy vertical scanning 
(FIGURE 8) will be performed. As soon as computer 112 first detects a characteristic change in 
the amplitude of penetrating ultrasonic pulses, it means that the pulses have first met mandibular 

15 canal 38 on their path from first transducer 106 to second transducer 108. Subsequently, when 
amplitude returns to its level prior to the characteristic change, it means that the penetrating 
ultrasonic pulses no longer pass through mandibular canal 38. Computer 112 then calculates and 
displays the diameter of mandibular canal 38 based on the detected borders of the canal and on 
location information supplied by location monitoring device 116. 

20 It is yet another particular feature of apparatus 100 that it can detect and locate an internal 

structure within hard tissue, for example mandibular canal 38 within mandible 30, based on 
spectral function analysis. Transducers 106 and 108 scan mandible 30, and ADC 110 sends to 
computer 112 samples of the received ultrasonic pulses, all as explained hereinabove. Computer 
112 processes the samples and produces, for each pair of locations of the transducers, a spectral 

25 function representing ultrasonic pulses emitted from and received at those locations. The spectral 
functions can be produced, for example, using a known Fast Fourier Transform (FFT) algorithm. 
The computer records, for each given location of the transducers, a set containing the spectral 
function of the received ultrasonic pulses along with the corresponding location of the 
transducers. Computer 112 can then analyze various characteristics of the spectral functions of 

30 the received ultrasonic pulses to determine the location of an internal structure within the hard 
tissue being examined. For example, it is known in the art that if an ultrasonic signal propagates 
through a medium which causes attenuation, the energy of the signal will be concentrated in a 
lower frequency range than in the original signal. In general, a higher level of attenuation will 
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cause a greater frequency shift. As mentioned heretofore, the present invention teaches that 
mandibular canal 38 often causes a different level of attenuation than cortical 34 and trabecular 
bone 36 around it. It therefore follows that the frequency spectrum of ultrasonic pulses received 
by second transducer 108 afler having traversed mandible 30 via mandibular canal 38 
5 (FIGURES 8c and 9b), will often be shifted compared to the frequency spectrum of similar 
ultrasonic pulses which did not travel through the mandibular canal. 

FIGURE 10 illustrates an example of three spectral functions of the same ultrasonic 
signal: a first spectral function 210 represents the ultrasonic signal as it is originally emitted into 
a jaw; a second spectral function 212 represents the same ultrasonic signal as received after 

10 having propagated mainly through the trabecular bone of the jaw (FIGURES 8a, 8b or 9a); and a 
third spectral function 214 representing the same ultrasonic signal as received after having 
traveled partly through trabecular bone and partly through the mandibular canal (FIGURES 8c or 
9b). In this example, the frequencies of peak amplitudes P212 and P214 of spectral functions 212 
and 214 (i.e. of the attenuated received signals) are lower than the frequency of peak amplitude 

15 P210 of spectral function 210 (i.e. of the originally emitted signal). Furthermore, the example 
illustrates that the frequency of peak amplitude P214 of ultrasonic pulses which traveled partly 
through mandibular canal is higher than the frequency of peak amplitude P212 of ultrasonic 
pulses which traveled mainly through trabecular bone. Hence, when computer 112 detects a 
characteristic frequency shift in the spectral function of received ultrasonic pulses, it will notify 

20 the user (e.g. by issuing an indication on display 114) that an internal structure lies on the 
ultrasonic travel path between first transducer 106 and second transducer 108. 

It is yet another particular feature of apparatus 100 that it can detect and locate an internal 
structure within hard tissue, based on artificially amplified characteristics of the received pulses, 
such characteristics being indicative of the internal structure. In some cases, the effect of the 

25 internal structure of interest on ultrasonic pulses traveling through it, may not be easily 
distinguishable from the effect of the surrounding area. This is the case, for example, when the 
level of attenuation inside mandibular canal 38 is nearly equal to the level of attenuation inside 
cortical 34 and trabecular bone 36. In such cases, computer 112 may perform ftirther 
manipulations on the samples which represent the received ultrasonic pulses, in order to amplify 

30 certain characteristics which enable detecting the internal structure. 

According to this method, after computer 112 obtains spectral functions of ultrasonic 
pulses received in various locations across mandible 30 (as explained hereinabove), the computer 
processes the recorded spectral ftinctions to determine a discrete representative value for each 
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spectral function. Such a representative value can be, for example, the minimal, maximal, 
average or root mean square (RMS) amplitude of each spectral function, or the total sum of all 
amplitudes in each spectral function, or any other predetermined representation criterion. As a 
result, the computer obtains and records a first array of representative values as a function of the 
5 locations of the transducers across mandible 30. Following, the computer analyzes the first array, 
in order to detect a distinct local variability in the representative values. In case such distinct 
local variability is found, the computer can conclude and notify the user (e.g. on display 114) 
that the internal structure of interest is situated in the locations corresponding to the locally 
varying representative values. If no such distinct local variability is detected, the computer will 

10 proceed and obtain a second array of representative values, which is based on a different 
representation criterion than that of the first array. If the first array is, for example, an array of 
the total sum of amplitudes, then the second array can be, for example, an array of the RMS of 
amplitudes. The computer analyzes the second array to find a distinct local variability which 
indicates the presence of the internal structure of interest. If again no such distinct local 

1 5 variability is detected, the computer will perform a predetermined mathematical manipulation on 
the first and\or second arrays. For example, the computer can divide the first array by the second 
array to obtain a third, synthetic array of representative values as a function of locations across 
the mandible. The computer analyzes the third, synthetic array to find a distinct local variability 
which indicates the presence of the internal structure. Often the synthetic array will be 

20 characterized by a distinct local variability, even though the first and second arrays which 
formed the third array were not so characterized. This is because certain mathematical 
manipulations may amplify physical effects which indicate the presence of the internal structure 
of interest. The computer may repeat the above procedure a desired number of repetitions. 

It is yet another particular feature of appeuatus 100 that it enables calculating the diameter 

25 of an internal structure within hard tissue, based on analysis of the ratio between the amplitude of 
received ultrasonic pulses and the attenuation coefficients inside the organ being examined. 
According to this method, after computer 112 obtains spectral functions of ultrasonic pulses 
received in various locations across mandible 30 (as explained hereinabove), the computer 
determines the minimal amplitude A| of ultrasonic pulses received by second transducer 108 

30 during the scanning process, and the maximal amplitude Ai of such pulses. 

As previously taught herein, attenuation in mandibular canal 38 is typically different than 
in cortical 34 and trabecular bone 36. Where attenuation in the mandibular canal is lower than in 
cortical and trabecular bone, minimal amplitude A| will be measured when penetrating ultrasonic 
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pulses do not travel through mandibular canal 38 (FIGURES 8a, 8b or 9a), whereas maximal 
amplitude A2 will be measured when the pulses travel partly via mandibular canal 38 (FIGURES 
8c or 9b). Horizontal width of trabecular area 36 is approximately equal to the total horizontal 
travel distance of ultrasonic pulses (in case of vertical scanning). Similarly, the vertical width Sy 
5 is approximately equal to the total vertical travel distance of ultrasonic pulses (in case of 
horizontal scanning). The attenuation coefficient inside trabecular area 36, designated 71, as well 
as the attenuation coefficient inside mandibular canal 38 (approximately equal to the attenuation 
coefficient in water) designated 72, are both measured using conventional means and are fed into 
computer 112. 

10 It is known in the art that the ratio between the amplitude A of a received ultrasonic 

signal and the distance S traveled by the ultrasonic signal inside a medium with an attenuation 
coefficient y can be expressed by the following first equation: 

A « e-<'-^> 

Hence, the ratio between minimal amplitude Ai and trabecular width S, can be expressed 
15 by the following second equation: 

And the ratio between maximal amplitude A2 and diameter D, of mandibular canal 38 
can be expressed by the following third equation: 

« ^-(ri (5x-^x)+(r2-^x)) 

20 If the third equation is divided by the second equation, then the following fourth equation 

is obtained: 

A 

The above fourth equation can be solved for Dx as expressed by the following fifth 
equation: 
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D.= ^ 

Thus, having acquired minimal and maximal amplitude measurements A| and A2 and 
attenuation coefficients yi and 72, computer 112 calculates and displays on display 114 horizontal 
diameter of mandibular canal 38, using the above fifth equation. Vertical diameter Dy can be 
5 measured in a similar manner, based on minimal and maximal amplitudes Ai and A2 measured 
vertically. 

It is yet another particular feature of apparatus 100 that it can detect and locate an internal 
structure within hard tissue (e.g. sinal and\or nasal cavities in a maxilla), based on analysis of the 
ratio between the travel distance and attenuation of ultrasonic pulses within the hard tissue. 

10 Referring now to FIGURE 15, this method will now be described. As explained hereinabove, 
transducers 106 and 108 are coupled using coupling material 74 to opposite surfaces of the organ 
being examined whilst facing each other. For example, in FIGURE 15 first transducer 106 is 
coupled to a first surface 52 of a maxilla 50, and second transducer 108 is coupled to a second 
surface 52' thereof In order to ensure that the transducers are facing each other (e.g. when the 

15 surface to which first transducer 106 is coupled is not parallel to the surface to which second 
transducer 108 is coupled), angular transducers may be used. When the transducers are properly 
coupled, the user starts the operation of apparatus 100 using controller 102. The controller 
commands location monitoring device 116 to determine the initial locations of the transducers, 
and to further determine the travel distance TD of ultrasonic pulses inside maxilla 50. Location 

20 monitoring device 116 determines the abovementioned factors using conventional measuring 
means and sends the measurements to computer 112. Controller 102 further commands 
pulser\receiver 104 to generate electrical pulses and to send them to first transducer 106 and to 
ADC 110. First transducer 106 converts the electrical pulses to ultrasonic pulses and emits them 
toward first surface 52. ADC 110 samples the electrical pulses which represent the emitted 

25 pulses and sends these samples to computer 112, The emitted ultrasonic pulses partially penetrate 
first surface 52 and then travel inside maxilla 50 eventually reaching second transducer 108. 
Second transducer 108 receives ultrasonic pulses which propagated through maxilla 50, converts 
them into analog electrical pulses and sends the electrical pulses to ADC 110 via pulserXreceiver 
104. The ADC samples the electrical pulses which represent the received pulses and sends the 
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samples to computer 112. The computer records, as a first set of results, the samples which 
represent the emitted pulses and the samples which represent the received pulses along with the 
initial locations of the transducers and the initial travel distance TD. 

The user scans maxilla 50 according to a desired scanning pattern, for example, by 
5 moving first transducer 106 along first surface 52 whilst also moving second transducer 108 
along second surface 52' in parallel to the first transducer. As mentioned hereinabove, apparatus 
100 can be equipped with automatic scanning means which perform this action instead of the 
user, either mechanically or electronically. During the scanning process, ADC 110 samples the 
ultrasonic pulses emitted and received at each new location of the transducers, and computer 112 

10 records sets of samples, along with their corresponding transducer locations and travel distance 
TD. The computer calculates, for each new location of the transducers, a ratio R between 
ultrasonic attenuation and travel distance. 

According to the teachings of the present invention, ratio R will remain approximately 
the same as long as penetrating ultrasonic pulses do not encounter an internal structure which 

15 causes higher or lower attenuation thsm its proximate surroundings within the hard tissue, on 
their travel path from first transducer 106 to second transducer 108. However, if penetrating 
ultrasonic pulses encounter such an internal structure, then ratio R will necessarily change. 
FIGURE 15 illustrates an example in which ultrasonic pulses traveling inside maxilla 50 
encounter sinal and\or nasal cavities 60 on their travel path from first transducer 106 to second 

20 transducer 108. Since sinal and\or nasal cavities 60 are sometimes filled with air, they may cause 
higher attenuation than cortical 54 and trabecular bone 56. Thus, if ultrasonic pulses traverse 
cavities 60 via travel path TPi and a travel distance TDi, then they will suffer higher attenuation 
on their way to second transducer 108 than similar pulses traveling an equivalent travel distance 
TDi via cortical and trabecular bone only. Alternatively, if the ultrasonic pulses bypass cavities 

25 60 via travel path TP2, they will still suffer higher attenuation because their travel distance will 
be longer than TDi. In both scenarios, computer 112 will detect the internal structure (cavities 
60), based on the change in ratio R. When computer 112 detects such a change in ratio R, it will 
notify the user (e.g. by issuing an indication on display 114) that an internal structure lies on the 
path between first transducer 106 and second transducer 108. 

30 It is yet another particular feature of apparatus 100 that it can measure the depth of an 

internal structure within hard tissue, such as the depth of the sinal and\or nasal cavities in a 
maxilla, based on changes in ratio R. According to this method, the transducers scan the maxilla 
as explained heretofore. As soon as computer 112. detects a characteristic change in ratio R, it 
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means that the penetrating ultrasonic pulses have first met cavities 60 on their path from first 
transducer 106 to second transducer 108. Since the computer receives from location monitoring 
device 116 the location of the transducers at any given moment, the computer can therefore 
calculate vertical depth Ly of cavities 60 as the vertical distance from alveolar ridge 64 of 
5 maxilla 50 to the location of the transducers which first introduced a characteristic change in 
ratio R. The measured depth will then be displayed on display 114. 

It is still another particular feature of apparatus 100 that it can combine some or all of the 
above methods to produce an integrated internal image of hard tissue. For example, computer 
112 in apparatus 100 can produce and display on display 114 a sectional image of mandible 30 

10 based on such parameters as vertical depth Ly, horizontal depth L,, vertical diameter Dy and 
horizontal diameter D, of mandibular canal 38 (all obtained using the teachings of the present 
invention), combined with such additional parameters as vertical width Sy and horizontal width 
Sx of mandible 30 (obtained using conventional measuring means). Subsequently, computer 112 
can further integrate several such sectional images of mandible 30 to create and display a 

1 5 three-dimensional internal image of the mandible. 

Referring now to FIGURE 11, an alternative embodiment of the present invention, 
designated 200, will be described. Apparatus 200 is an apparatus for non-invasive ultrasonic 
imaging of hard tissue. Apparatus 200 is partially similar to apparatus 100, and therefore 

20 common elements will be hereunder denoted with the same reference numerals. Apparatus 200 
comprises a controller 102 coupled to a pulser\receiver 104 and to a location monitoring device 
116. PulserXreceiver 104 is coupled to a transducer 120 capable of emitting and receiving 
ultrasound, and to an analog-to-digital converter (ADC) 110. ADC 110 is coupled to a computer 
112, which is coupled to a display device 114. Location monitoring device 116 is coupled to 

25 transducer 120 and to computer 112. 

Referring now to FIGURE 12 the method of operation of apparatus 200 will be 
described. Apparatus 200 operates according to a pulse-echo method, as opposed to apparatus 
100 which is based on a through-transmission method. A user couples transducer 120 to an initial 
location on a surface of an organ being examined, for example, on first surface 32 of mandible 

30 30. In FIGURE 12, first surface 32 is, for example, a toothless buccal gum surface. As explained 
hereinabove in reference to transducers 106 and 108 of apparatus 100, the same procedure can be 
performed on other suitable surfaces inside or outside the mouth cavity. In order to improve 
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penetration of ultrasound signals into the jaw, the user applies coupling material 74 between 
transducer 120 and surface 32. 

When transducer 120 is properly coupled, the user starts the operation of apparatus 200 
using controller 102. The controller commands location monitoring device 116 to determine the 
5 initial location of transducer 120. Device 116 detemiines the initial location using a conventional 
measuring method and sends the result to computer 112. The controller further commands 
pulserVeceiver 104 to generate electrical pulses and to send them to transducer 120. Transducer 
120 converts the electrical pulses to ultrasonic pulses and emits them toward first surface 32. The 
emitted ultrasonic pulses partially penetrate first surface 32 via coupling material 74. Some of 

10 the penetrating ultrasonic pulses then travel inside mandible 30 through cortical bone layer 34 
and trabecular bone 36. Most of the penetrating pulses travel inside mandible 30 in an 
approximately straight trajectory (depicted in the drawings as a broken-line arrow). The 
penetrating ultrasonic pulses are then partially reflected back toward transducer 120 from the 
opposite cortical bone layer 34*. In some cases, the penetrating pulses may also travel through 

15 mandibular canal 38 on their way to opposite surface 34' and after being reflected back toward 
transducer 120. Transducer 120 receives ultrasonic echoes reflected from within mandible 30, 
converts them into analog electrical pulses and sends the electrical pulses to ADC 110 via 
pulser\receiver 104. The ADC samples the electrical pulses which represent the received 
ultrasonic echo pulses and sends these samples to computer 112. The computer records, as a first 

20 set of results, the samples which represent the received pulses along vwth the initial location of 
the transducer. 

The user scans mandible 30, for example, by moving transducer 120 vertically and\or 
horizontally along the surface of mandible 30. FIGURES 13a - 13c illustrate an example of a 
vertical scanning movement- FIGURE 13a shows the initial location of transducer 120 which is 

25 relatively far from mandibular canal 38. FIGURE 13b shows the location of the transducer after 
having been moved approximately halfway closer to mandibular canal 38. FIGURE 1 3c shows 
the location of the transducer after having been moved fiirther down, so that mandibular canal 38 
lies on the ultrasonic travel path (depicted as a broken-line arrow) extending straightforward 
from the transducer. FIGURES 14a - 14b illustrate an example of a horizontal scanning 

30 movement of transducer 120 along alveolar ridge 40. FIGURE 14a shows the initial location of 
transducer 120. The transducer is initially placed so that the travel path (depicted as a broken-line 
arrow) of ultrasonic pulses and echoes does not pass through mandibular canal 38. FIGURE 14b 
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shows the transducer after having been moved toward the middle of the jaw, so that mandibular 
canal 38 lies on the ultrasonic travel path extending straightforward from the transducer. 

Alternatively, the user may scan mandible 30 by moving the transducer in a sweeping 
motion, or by any other desired scanning pattern. Apparatus 200 can further be equipped with 
5 automatic mechanical scanning means, for example, a slidable mount to which transducer 120 is 
attached and which is programmed to execute a desired scanning pattern. Alternatively, 
apparatus 200 can be equipped with electronic scanning means, for example, an array of 
transducers instead of single transducer 120. 

During the scanning process, ADC 110 samples the ultrasonic pulses received at each 

10 new location of the transducers, and computer 112 keeps recording sets of samples along with 
the corresponding location of the transducer. Computer 112 is programmed to detennine certain 
physical characteristics of the hard tissue being examined, based on an analysis of the recorded 
sets. As mentioned hereinabove in reference to apparatus 100, the present invention teaches that 
when, during the scanning process, mandibular canal 38 is situated on the path extending 

15 straightforward from transducer 120, consequently computer 112 will detect a local change 
(increase or decrease) in amplitude of the penetrating pulses. 

Hence, it is a particular feature of apparatus 200 that it detects and locates an internal 
structure within hard tissue (e.g. the mandibular ceuial in a mandible), by detecting a local change 
in the amplitude of ultrasonic echo pulses. When computer 112 detects such a local change in the 

20 amplitude, it will notify the user (e.g. by issuing an indication on display 114) that the internal 
structure lies on the path extending straightforward from transducer 120. 

It is another particular feature of apparatus 200 that it can measure the depth of an 
internal structure within hard tissue, for example, the depth of the mandibular canal in a 
mandible, based on amplitude measurements. Transducer 120 scans mandible 30 according to 

25 the method described above. As soon as computer 112 first detects a characteristic change in the 
amplitude of received ultrasonic echo pulses, it means tliat the pulses have first met mandibular 
canal 38 on their path from first transducer 120 to the surface of reflection. Since the computer 
receives from location monitoring device 116 the location of the transducer at any given 
moment, the computer can therefore calculate the depth of mandibular canal as follows. In case 

30 of vertical scanning (FIGURES 13a - 1 3c), vertical depth Ly will be calculated as the vertical 
distance from alveolar ridge 40 to the location of the transducer which first introduced a 
characteristic change in amplitude. Similarly, in case of horizontal scanning (FIGURES 14a - 
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14b), horizontal depth will be calculated as the horizontal distance from surface 32 to the 
location of the first characteristic amplitude change. 

It is yet another particular feature of apparatus 200 that it can measure the diameter of an 
internal structure of interest within hard tissue, for example, the diameter of mandibular canal 38 
5 within mandible 30, based on amplitude measurements. Transducer 120 scans mandible 30 
according to the method described hereinabove. When measuring horizontal diameter Dx 
horizontal scanning (FIGURE 14) will be performed, and when measuring vertical diameter Dy 
vertical scanning (FIGURE 13) will be performed. As soon as computer 112 first detects a 
characteristic change in the amplitude of received ultrasonic echo pulses, it means that the pulses 

10 have first met mandibular canal 38 on their path from transducer 120 to the surface of reflection. 
Subsequently, when amplitude returns to its level prior to the characteristic change, it means that 
the ultrasonic pulses no longer pass through mandibular canal 38. Computer 112 then calculates 
and displays the diameter of mandibular canal 38 based on the detected borders of the canal and 
on location information supplied by location monitoring device 116. 

15 It is yet another particular feature of apparatus 200 that it can detect and locate an internal 

structure within hard tissue, for example mandibular canal 38 within mandible 30, based on 
spectral function analysis. Transducer 120 scans mandible 30, and ADC 110 sends to computer 
112 samples of the received ultrasonic echo pulses, all as explained hereinabove. Computer 112 
processes the samples and produces (e.g. using a known PPT algorithm), for each given location 

20 of the transducer, a spectral function representing ultrasonic pulses emitted and received at that 
location. The computer records, for each given location of the transducer, a set containing the 
spectral function of the received ultrasonic pulses along with the corresponding location of the 
transducer. Computer 112 can then analyze various characteristics of the spectral functions of the 
received ultrasonic pulses to determine the location of an internal structure wdthin the hard tissue 

25 being examined. As mentioned heretofore in reference to apparatus 100, a frequency shift in the 
frequency spectrum of the received ultrasonic pulses may indicate the presence of mandibular 
canal 38 within mandible 30. Hence, when computer 112 detects a characteristic frequency shift 
in the spectral function of received ultrasonic echo pulses, it will notify the user (e.g. by issuing 
an indication on display 114) that an internal structure lies on the ultrasonic travel path extending 

30 straightforward from transducer 120. After apparatus 200 detects and locates the internal 
structure of interest, the apparatus may further be used in order to measure the depth of such 
internal structure using location information supplied by location monitoring device 116, as 
explained heretofore. 
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It is yet another particular feature of apparatus 200 that it can detect and locate an internal 
structure within hard tissue, based on artificially amplified characteristics of the received pulses, 
such characteristics being indicative of the internal structure. After computer 112 obtains spectral 
functions of ultrasonic pulses received in various locations along mandible 30 (as explained 
5 hereinabove), the computer processes the recorded spectral functions to determine a discrete 
representative value for each spectral function. As a result, the computer obtains and records a 
first array of representative values as a function of the locations of the transducer across 
mandible 30. Following, tlie computer analyzes the first array, in order to detect a distinct local 
variability in the representative values. In case such distinct local variability is found, the 

10 computer can conclude and notify the user (e.g. on display 114) that the internal structure of 
interest is situated in the locations corresponding to the locally varying representative values. If 
no such distinct local variability is detected, the computer will proceed and obtain a second array 
of representative values, which is based on a different representation criterion than that of the 
first array. If the first array is, for example, an array of the total sum of amplitudes, then the 

1 5 second array can be, for example, an array of the RMS of amplitudes. The computer analyzes the 
second array to find a distinct local variability which indicates the presence of the internal 
structure of interest. If again no such distinct local variability is detected, the computer will 
perform a predetermined mathematical manipulation on the first and\or second arrays. For 
example, the computer can divide the first array by the second array to obtain a third, synthetic 

20 array of representative values as a function of locations across the mandible. The computer 
analyzes the third, synthetic array to find a distinct local variability which indicates the presence 
of the internal structure. The computer may repeat the above procedure a desired number of 
repetitions. 

It is yet another particular feature of apparatus 200 that it enables calculating the diameter 
25 of an internal structure within hard tissue, based on analysis of the ratio between the amplitude of 
received ultrasonic pulses and the attenuation coefficients inside the organ being examined. 
According to this method, after computer 112 obtains spectral functions of ultrasonic pulses 
received in various locations along mandible 30 (as explained hereinabove), the computer 
determines the minimal amplitude A3 of ultrasonic pulses received by transducer 120 during the 
30 scanning process, and the maximal amplitude A4 of such pulses. As mentioned hereinabove, the 
present invention teaches that attenuation in mandibular canal 38 is typically different than in 
cortical 34 and trabecular bone 36. Where attenuation in the mandibular canal is lower than in 
cortical and trabecular bone, minimal amplitude A3 will be measured when penetrating ultrasonic 
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pulses do not travel through mandibular canal 38 (FIGURES 13a, 13b or 14a), whereas maximal 
amplitude A4 will be measured when the pulses travel partly via mandibular canal 38 (FIGURES 
13c or I4b). The attenuation coefficient inside trabecular bone 36, designated yi, as well as the 
attenuation coefficient inside mandibular canal 38 (approximately equal to the attenuation 
coefficient in water) designated Y2» are both measured using conventional means and are fed into 
computer 112. 

As mentioned above, it is known in the art that the ratio between the amplitude A of an 
ultrasonic signal and the distance S traveled by the ultrasonic signal inside a medium with an 
attenuation coefficient y can be expressed by the following first equation: 



Hence, the ratio between minimal amplitude A3 and trabecular width S, can be expressed 
by the following sixth equation: 

The ratio between maximal amplitude A4 and diameter Dx of mandibular canal 38 can be 
15 expressed by the following seventh equation: 



If the seventh equation is divided by the sixth equation^ then the following eighth 
equation is obtained: 

A _ g2 Z), .(r,-r2) 

A 

20 The above eighth equation can be solved for Dx as expressed by the following ninth 

equation: 

A. 



log 



X 



Thus, having acquired minimal and maximal amplitude measurements Aj and A4, and 
attenuation coefficients yi and y2, computer 112 calculates and displays on display 114 horizontal 
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diameter Dx of mandibular canal 38, using the above ninth equation. Vertical diameter Dy can be 
measured in a similar manner, based on minimal and maximal amplitudes A3 and A4 measured 
vertically. 

It is still another particular feature of apparatus 200 that it can combine some or all of the 
5 above methods to produce an integrated internal image of hard tissue. For example, computer 
112 in apparatus 200 can produce and display on display 114 a sectional image of mandible 30, 
based on such parameters as vertical depth Ly. horizontal depth Lx, vertical diameter Dy and 
horizontal diameter Dx of mandibular canal 38 (all obtained using the teachings of the present 
invention), combined with such additional parameters as vertical width Sy and horizontal width 
10 Sx of mandible 30 (obtained using conventional measuring means). Subsequently, computer 112 
can further integrate several such sectional images of mandible 30 to create and display a 
three-dimensional internal image of the mandible. 

Referring now to FIGURE 16, yet another alternative embodiment of the present 

15 invention, designated 300, will be described. Apparatus 300 is an apparatus for non-invasive 
ultrasonic imaging of hard tissue. Apparatus 300 is partially similar to apparatuses 100 and 200, 
and therefore common elements will be hereunder denoted with the same reference numerals. 
Apparatus 300 comprises a controller 102 coupled to a pulserXreceiver 104 and to a location 
monitoring device 116. PulserVreceiver 104 is coupled to a multiplexer (MUX) 130, which is 

20 coupled to a first transducer 106 and to a second transducer 108. Transducers 106 and 108 are 
capable of emitting and receiving ultrasound. PulserVeceiver 104 is fxirther coupled to an 
analog-to-digital converter (ADC) 110, coupled to a computer 112, which is coupled to a display 
device 114. Location monitoring device 116 is coupled to first transducer 106, to second 
transducer 108, and to computer 112, 

25 Apparatus 300 is capable of operating in several modes of operation. In a first mode of 

operation, the through-transmission mode, a user sets MUX 130 so that it directs an output 
electrical signal from pulser\receiver 104 to first transducer 106, and further directs the analog 
electrical signal which represents the ultrasonic signal received by second transducer 108 back to 
the pulserXreceiver 104. In this mode apparatus 300 actually operates like apparatus 100, as 

30 explained hereinabove. In a second mode of operation, the pulse-echo mode, the user sets MUX 
130 so that it directs an output electrical signal from pulserXrecciver 104 to first transducer 106, 
and further directs the analog electrical signal which represents the ultrasonic signal received by 
first transducer 106 back to the pulserXreceiver. In this mode apparatus 300 actually operates like 
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apparatus 200, as explained hereinabove. In a third mode of operation, a combined mode, 
apparatus 300 operates in the first mode to obtain a first set of measurements; and then apparatus 
300 operates in the second mode to obtain a second set of measurements; computer 112 records 
the first and second sets of measurements, and integrates them to produce and display on display 
5 114 an integrated internal image of the hard tissue being examined, for example, a sectional 
image or a three-dimensional image based on a plurality of sectional images. 

Referring now to FIGURE 17 still another alternative embodiment of the present 
invention, designated 400, will be described. Apparatus 400 is an apparatus for non-invasive 

10 ultrasonic imaging of hard tissue. Apparatus 400 is partially similar to apparatus 300, and 
therefore common elements will be hereunder denoted with the same reference numerals. 
Apparatus 400 comprises a controller 102 coupled to a pulserVreceiver 104 and to a location 
monitoring device 116. PulserVreceiver 104 is coupled to a multiplexer (MUX) 130 and also to an 
analog-to-digital converter (ADC) 110. ADC 110 is coupled to a computer 112 coupled to a 

15 display device 114. MUX 130 is further coupled to an array of transducers 140 containing a 
plurality of transducers, for example, eight transducers 141 - 148. Each transducer in array 140 
is capable of emitting and receiving ultrasound. Apparatus 400 further comprises a location 
monitoring device 116 coupled to each of the transducers in array 140 and further coupled to 
computer 112. 

20 Referring now to FIGURE 18, the method of operation of apparatus 400 will be 

described. A user couples transducers 141 - 148 to an organ being examined so that at least some 
of the transducers are facing each other. For example, transducers 141 — 148 may be placed 
around a gum surface 32 of a mandible 30 as illustrated in FIGURE 18. The user applies a 
coupling material 74 between each transducer and the surface to which it is coupled, as explained 

25 hereinabove in reference to apparatus 100. When the transducers are properly coupled, the user 
starts the operation of apparatus 400 using controller 102. The controller commands location 
monitoring device 116 to determine the position of each of the transducers in array 140 with 
respect to the other transducers therein and with respect to surface 32 to which they are coupled. 
Device 116 sends the locations to computer 112 which records them, 

30 Next, apparatus 400 electronically scans the organ being examined, as follows. Controller 

102 commands pulser\receiver 104 to generate electrical pulses and send them to MUX 130 and 
to ADC 110. The ADC samples the electrical pulses, and sends to computer 112 the samples 
which represent the emitted pulses. MUX 130 directs the electrical pulses from pulser\receiver 



wo 02/085178 



PCT/IL02/00311 



30 

104 to a first transducer 141. The ultrasonic pulses emitted from first transducer 141 partially 
penetrate mandible 30 and then travel inside the mandible eventually reaching one or more of the 
remaining transducers in array 140. The receiving transducers convert the received ultrasonic 
pulses to analog electrical pulses, and send the latter to MUX 130. The MUX directs analog 
5 electrical pulses from the remaining transducers back to ADC 110 via pulserVreceiver 104. The 
ADC samples the electrical pulses and sends the samples, which represent the received ultrasonic 
pulses, to computer 112. The computer records the digital samples in discrete sets, wherein each 
set contains samples representing ultrasonic pulses emitted by one certain transducer and 
received by another certain transducer, along with the locations of these transducers (as 

10 determined by device 116). The above procedure is repeated several times, whilst each time a 
different transducer from array 140 emits ultrasound and the remaining transducers receive 
ultrasound. FIGURE 18 illustrates an example of various possible ultrasonic travel paths 
(depicted as broken lines) inside mandible 30. 

At the end of the scanning process described heretofore, computer 112 will have 

15 accumulated a desired number of sets of samples and travel paths. The computer processes the 
accumulated sets, and converts each set of samples and ultrasonic travel path into a set of a 
predetermined integral physical quantity and travel path. For example, the computer can 
calculate, based on samples of the emitted and received ultrasonic pulses in each travel path, the 
integral attenuation in each travel path. Following, the computer implements a known method, 

20 such as a Radon Transform algorithm, in order to deduce from the sets of integral physical 
quantity and travel path, a physical quantity per each point inside the hard tissue. For example, 
the computer can use a Radon Transform algorithm for deducing from all the sets of integral 
attenuation per travel path, the level of attenuation in each discrete point inside the scanned 
cross-section of the mandible. The computer can then display on display 114 a sectional internal 

25 image of mandible 30, based on the deduction algorithm. 

It is a particular feature of apparatus 400 that it can detect, locate and measure the size of 
an internal structure within hard tissue, for example, a mandibular canal 38 in a mandible 30. As 
mentioned above, the level of ultrasonic attenuation inside mandibular canal 38 is typically 
different than in cortical 34 and trabecular bone 36. The level of attenuation at each point inside 

30 mandible 30 can be determined using apparatus 400 as explained hereinabove. Following, 
computer 112 of apparatus 400 can detect mandibular canal 38 by finding the locations within 
mandible 30 that are characterized with an attenuation level other than that of its surroundings. 
Horizontal depth Lx, vertical depth Ly, horizontal diameter D,, and vertical diameter Dy (see 
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FIGURE 12) of mandibular canal 38, as well as other parameters pertaining to the size and 
location of the internal structure of interest, can all be calculated based on location information 
obtained by device 116 and on the Radon Transform internal image of mandible 30. * 

It is another particular feature of apparatus 400 that it can further integrate several 
5 sectional images of the hard tissue being examined, to create and display a tliree-dimensional 
internal image of the hard tissue. 

Thus, it is evident that the present invention provides a real-time, chair-side, accurate, 
safe, radiation-free, and economical method and apparatus for non-invasive ultrasonic imaging 
of hard tissue. 

10 While preferred embodiments of the present invention have been disclosed hereinabove, 

it is to be understood that these preferred embodiments are given as an example only and are not 
intended to be limiting. Those skilled in the art may make various modifications and additions to 
the embodiments used to illustrate the teachings of the present invention and those modifications 
and additions would remain within the scope of the present invention. 

15 It is emphasized that the present invention is not limited to imaging human jaws. Other 

uses, including but not limited to imaging other types of human hard tissue, in addition to animal 
hard tissue, are included in the scope of the present invention. 

Although the invention has been described in conjunction with specific embodiments 
thereof, it is evident that many alternatives, modifications and variations will be apparent to 

20 those skilled in the art. Accordingly, it is intended to embrace all such alternatives, 
modifications and variations that fall within the spirit and broad scope of the appended claims. 



wo 02/085178 



PCT/IL02/00311 



10 



15 



32 

WHAT IS CLAIMED IS: 

1 . A method of ultrasonic imaging of a biological tissue, comprising the steps of: 

(a) scanning the biological tissue with ultrasonic energy transmitted from and 
received at a plurality of known transmittal and reception locations about the 
biological tissue, to obtain a corresponding plurality of sets of digital data; and 

(b) processing said sets of digital data to produce an image of a cross-section of the 
biological tissue. 

2. The method of claim 1, wherein the biological tissue includes a jaw, 

3. The method of claim 1 , wherein the biological tissue includes a bone. 

4. The method of claim 3, wherein said bone is selected from the group consisting of a 
mandible and a maxilla. 

5. The method of claim 1, wherein said step of scanning includes applying at least one 
scanning mode selected from the group consisting of through-transmission mode and 
pulse-echo mode. 

20 6. The method of claim 1, wherein said step of scanning includes using at least one scanning 

mechanism selected from the group consisting of mechanical scanning and electronic 
scanning 

7. The method of claim 1 , wherein said step of scanning includes the substeps of: 
25 (i) disposing at least one ultrasonic transducer in proximity to the biological tissue; 

(ii) determining a location of said at least one ultrasonic transducer; 

(iii) transmitting ultrasonic energy into the biological tissue using said at least one 
ultrasonic transducer, said ultrasonic energy being reflected from the biological 
tissue; and 

30 (iv) receiving said ultrasonic energy reflected from the biological tissue using said at 

least one ultrasonic transducer. 



8. The method of claim 7, further including the substeps of: 
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(v) changing said location of said at least one ultrasonic transducer; and 

(vi) repeating steps (ii) through (v) for a desired number of repetitions. 

9. The method of claim 1 wherein said step of scanning includes the substeps of: 

(i) disposing a plurality of ultrasonic transducers in proximity to the biological 
tissue; 

(ii) determining a location of each of said plurality of ultrasonic transducers; 

(iii) transmitting ultrasonic energy into the biological tissue, using at least one first 
ultrasonic transducer selected from said plurality of ultrasonic transducers, said 
ultrasonic energy propagating through the biological tissue; and 

(iv) receiving said through propagating ultrasonic energy, using at least one second 
ultrasonic transducer selected from said plurality of ultrasonic transducers, said 
at least one second ultrasonic transducer being different from said at least one 
first ultrasonic transducer. 



1 0. The method of claim 9, further including the substeps of: 

(v) varying the selection of at least one ultrasonic transducer from the group 
consisting of said at least one first ultrasonic transducer and said at least one 
second ultrasonic transducer; and 

(vi) repeating steps (iii) through (v) for a desired number of repetitions. 

1 1 . The method of claim 1, wherein said step of scanning includes the substeps of: 

(i) disposing a plurality of ultrasonic transducers in proximity to the biological 
tissue; 

(ii) determining a location of each of said plurality of ultrasonic transducers; 

(iii) transmitting ultrasonic energy into the biological tissue, using at least one 
ultrasonic transducer selected from said plurality of ultrasonic transducers, said 
ultrasonic energy being reflected from the biological tissue; 

(iv) receiving said ultrasonic energy reflected from the biological tissue using semie 
said at least one selected ultrasonic transducer used for transmitting; 

(v) varying the selection of said at least one selected ultrasonic transducer used for 
transmitting and receiving said ultrasonic energy; and 

(vi) repeating steps (iii) through (v) for a desired number of repetitions. 
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12. The method of claim 1, wherein said step of processing includes deriving, from each of 
said sets of digital data, at least one value related to an ultrasonic property of the 
biological tissue. 

13. The method of claim 12, wherein said ultrasonic property includes ultrasonic attenuation. 

14. The method of claim 1, wherein said step of processing includes the substep of: 

(i) generating, for each of said sets of digital data, a spectral function of the 
corresponding said received ultrasonic energy. 

15. The method of claim 14, wherein said substep of generating a spectral function includes 
applying a fast Fourier transform (FFT) algorithm. 

16. The method of claim 14, further including the substep of: 

(ii) representing each said spectral function by a first representative value, all such 
said first representative values forming a first array of representative values. 

17. The method of claim 16, wherein said first representative vialue is selected from the group 
consisting of a minimum, a maximum, an average, a root mean square (RMS), and a total 
sum of amplitudes. 

18. The method of claim 16, further including the substeps of: 

(iii) identifying at least one distinct local variability in said first array of 
representative values; and 

(iv) calculating the value of a dimension related to an internal structure within the 
biological tissue, based on said at least one distinct local variability. 

19. The method of claim 18, wherein said dimension related to an internal structure within 
the biological tissue is selected from the group consisting of a width, a length, a depth, a 
height, a diameter and a thickness. 

20. The method of claim 1 8, wherein said internal structure includes a cavity. 
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21. The method of claim 20, wherein said cavity is selected from the group consisting of a 
mandibular canal, a sinal cavity and a nasal cavity. 

5 22. The method of claim 16, further including the substeps of: 

(iii) representing each said spectral function by a second representative value, all 
such said second representative values forming a second array of representative 
value; and 

(iv) producing a synthetic array of representative values, based on said first and 
10 second arrays of representative values. 

23. The method of claim 22, further including the substeps of: 

(v) identifying at least one distinct local variability in said synthetic array of 
representative values; and 

15 (vi) calculating the value of a dimension related to an internal structure within the 

biological tissue, based on said at least one distinct local variability. 

24. The method of claim 23, wherein said dimension related to an internal structure within 
the biological tissue is selected from the group consisting of a width, a length, a depth, a 

20 height, a diameter and a thickness. 

25. The method of claim 23, wherein said internal structure includes a cavity. 

26. The method of claim 25, wherein said cavity is selected from the group consisting of a 
25 mandibular canal, a sinal cavity and a nasal cavity. 

27. The method of claim 1, wherein said step of processing includes applying a Radon 
transform algorithm. 

30 28, The method of claim 1 6, further including the substeps of: 

(iii) identifying at least one extreme value in said first array of representative values; 

(iv) measuring a first attenuation coefficient inside the biological tissue; 
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(v) measuring a second attenuation coefTicient inside an internal structure within 
the biological tissue; 

(vi) measuring the value of a first dimension related to the biological tissue; and 

(vii) calculating the value of a second dimension related to said internal structure, 
based on said at least one identified extreme value, said first attenuation 
coefficient, said second attenuation coefficient, and said value of said first 
dimension. 

29.. The method of claim 1, wherein said step of processing includes the substeps of: 

(i) deriving, from each of said sets of digital data, a property value related to an 
ultrasonic property of the biological tissue; 

(ii) measuring, for each of said sets of digital data, a dimension value of a 
dimension related to the biological tissue and to same said set of digital data; 
and 

(iii) calculating, for each of said sets of digital data, a set ratio between said property 
value and said dimension value. 

30. The method of claim 29, wherein said ultrasonic property includes ultrasonic attenuation. 

31. The method of claim 29, wherein said dimension related to the biological tissue is 
selected from the group consisting of a width, a length, a depth, a height, a diameter and a 
thickness. 

32. The method of claim 29, wherein said step of processing further includes the substeps of: 

(iv) forming an array of ratios based on said set ratios calculated for all said sets of 
digital data; 

(v) identifying at least one distinct local variability in said array of ratios; and 

(vi) calculating the value of a dimension related to an internal structure within the 
biological tissue, based on said at least one identified distinct local variability. 



33. The method of claim 32, wherein said internal structure includes a cavity. 
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34. The method of claim 33, wherein said cavity is selected from the group consisting of a 
mandibular canal, a sinal cavity and a nasal cavity. 

35. The method ofclaim 1, further comprising the steps of: 

5 (c) repeating steps (a) and (b) for a desired number of repetitions, to accumulate 

various cross-sectional images of the biological tissue; and 
(d) transforming said accumulated cross-sectional images into a three-dimensional 
image of the biological tissue. 

10 36. An apparatus for ultrasonic imaging of a biological tissue, comprising: 

(a) means for scanning the biological tissue with ultrasonic energy from a plurality of 
transmittal and reception locations about the biological tissue; 

(b) means for determining said transmittal and reception locations; 

(c) means for generating a plurality of sets of digital data, said sets of digital data 
15 representing parameters of said ultrasonic energy at each of said plurality of 

transmittal and reception locations; and 

(d) means for processing said plurality of sets of digital data to produce a cross- 
sectional image of the biological tissue. 

20 37. The apparatus of claim 36, wherein the biological tissue is selected from the group 

consisting of a jaw and a bone. 

38. The apparatus of claim 37, wherein said bone is selected from the group consisting of a 
mandible and a maxilla. 



25 



39. The apparatus of claim 36, wherein said means for scanning includes at least one means 
from the group consisting of means for scanning in through-transmission mode and 
means for scanning in pulse-echo mode. 



30 



40. The apparatus of claim 36, wherein said means for scanning includes: 

(i) a plurality of ultrasonic transducers for transmitting and receiving ultrasonic 
energy; and 
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(ii) means for changing the location of at least one of said plurality of ultrasonic 
transducers. 

41. The apparatus of claim 36, wherein said means for scanning includes: 

(i) a plurality of ultrasonic transducers for transmitting and receiving ultrasonic 
energy; and 

(ii) means for selecting at least one of said plurality of ultrasonic transducers. 
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FIG. 5a 
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FIG. 6 



100 



106 



FIRST 
TRANSDUCER 



108 



SECOND 
TRANSDUCER 



116 



TV 



LOCATION 
MONITORING 
DEVICE 



102 



CONTROLLER 



104 



PULSERX RECEIVER 



110 



A/D CONVERTER 



Ypl12 



COMPUTER 



_J___yp114 

DISPLAY 



wo 02/085178 



7/18 



PCT/IL02/00311 




wo 02/085178 



PCT/IL02/00311 



8/18 




wo 02/085178 



PCT/IL02/00311 



9/18 





wo 02/085178 



10/18 



PCT/IL02/00311 




wo 02/085178 



PCT/IL02/00311 



11/18 



FIG. 11 
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